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UNSTEADY OBLIQUE INTERACTION OF A SHOCK WAVE WITH A PLANE DISTURBANCE ^ 

By Franklin K. Moore 


SUMMARY 

Analysis is made of the flow field produced by oblique im- 
pingement of weak plane disturbances of arbitrary profile on a 
plane normal shock. Three types of disturbance are considered: 

{a) Sound wave propagating in the gas at rest into which the 
shock moves. The sound wave refracts either as a simple isen- 
tropic sound wave or as an attenuating isentropic pressure wave, 
depending on the angle between the shock and the incident sound 
wave. A stationary vorticity wave of constant pressure appears 
behind the shock. 

ib) Sound wave overtaking the shock from behind. The 
sound wave refiects as a sound wave, and a stationary vorticity 
wave is produced. 

(c) An incompressible vorticity wave stationary in the gas 
ahead of the shock. The incident wave refracts as a stationary 
vorticity wave, and either a sound wave or attenuating pressure 
wave is also produced. 

Computations are presented for the first two types of incident 
wave, over the range of incidence angles, for shock Mach num- 
bers of 1, 1.5, and » . 

INTRODUCTION 

The unsteady one-dimensional interaction of normal shock 
waves and disturbances, such as sound waves or other shock 
waves, has been studied quite thoroughly (an example is 
Kantrowitz’ paper on shock stabilit 3 r^ ref. 1). The steady 
interaction between normal shock waves and plane Mach 
waves has been treated by Adams (ref. 2). 

The general class of unsteady flow problems is currently 
of increasing interest, in connection particularly with sta- 
bility of high-speed aerodynamic and combustion processes. 
The effect of a shock passing through a flow field (or vice 
versa) is liliely to be important in many applications. For ex- 
ample, a hot-wire anemometer intended to measure the fluctuat- 
ingfield of turbulence in a supersonic stream will actuallymeas- 
ure the turbulence as modified by passage through the nearly 
steady bow shock of the probe. 

Considering, for simplicity, that the flow interacting with 
a normal shock is a nonviscous field of weak disturbance, it 
may usually be considered irrotational and isentropic (such 
as produced by a moving slender body) and therefore can be 
imagined to be composed of a suitable ai-ray of sound waves. 
Another possible type of weak nonviscous disturbance would 
be a stationary, incompressible flow of variable vorticity 


(turbulence which is convected rapidly past the point of 
observation is commonly thought of in this way). 

Either of these two types of flow may be represented as a 
linear combination of plane waves (each wave either a sound 
wave or a rotational wave, depending on the type of flow 
to be represented) of various amplitudes, wave lengths, 
and orientations. Thus, if the interactive effect of a shock 
and each constituent wave may be found by a linear analysis, 
the complete problem may in principle be solved by linear 
combination of the resulting flow fields behind the shock. 
The interaction between a turbulent field and a wind-tunnel 
screen or contraction, or both, has been successfully carried 
out in references 3 through 5 by this method. 

The present report concerns the interaction of a normal 
shock met obliquely by a plane sound wave or by a con- 
vected plane vorticity wave. Since sound waves may im- 
pinge on a shock either from upstream or downstream, both 
cases are considered. The oblique interaction of a shock 
and weak vorticity wave is also treated in a current investi- 
gation by Ribner (ref. 6). 

The shock is considered to be moving freely into gas nominally 
at rest (as in a shock tube, when wall effects are neglected) . 
Of course, if the observer moves at a constant speed with the 
shock, the flow appears as that associated with a steady 
shock, under dift'erent stagnation conditions. The shock- 
tube pointofviewis adopted in order that there be no question 
of how the equivalent steady shock is “anchored”; that is, 
end effects on the shock are not contemplated. 

GOVERNING UNSTEADY EQUATIONS 

In the following paragraphs, the equations will be derived 
which pertain to the propagation of a plane normal shock 
wave through a gas at rest, as modified by the influence of a 
weak pattern of unsteady disturbance. 

If the shock propagates without disturbance, its instan- 
taneous position is Xi=Vt, in a coordinate system fixed in 
the fluid nominally at rest ahead of the shock (fig. IJ. (AU 
symbols are defined in appendix A.) The constant velocity 
of the shock front is V, and the corresponding constant 
velocity of the gas behind the shock is U, in this system. 
This one-dimensional motion is considered to be perturbed 
slightly by the presence of a weak field of unsteady plane 
flow. The velocity ahead of the shock is written as Ui(x,y,t), 
Vi(x,y,t); behind the shock as U-[-Ui{x,y,f), V 2 {x,y,t). Pres- 
sure, density, and temperature are written as P-{-p, R f-p, 


* Supersedes NAO A TN 2879, “Unsteady Oblique Interaction of a Shock Wave with a Plane Disturbance,** by Franklin K. Moore, 1953. 
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and 04-0, respectively, where the capitalized symbols refer 
to the basic steady-shock motion, and the lower case to the 
unsteady disturbance. Throughout, subscripts 1 and 2 are 
used to specify conditions ahead of and behind the shock, 
respectively. As a result of the unsteady disturbance, the 
shock front itself undergoes a small unsteady displacement, 
given by ^{y,t). Thus, the position of the shock at any 
instant is 

SHOCK RELATIONS 

Because of the rapidity vuth which changes occur across 
a shock wave, the disturbed shock may be regarded as be- 
having in a locally quasi-steady manner; that is, in a coordi- 
nate system fixed in the shock at each instant, the usual 
Rankine-Hugoniot relations apply. Because of the disturb- 
ance, the shock is slightly oblique in such a coordinate system 
(see sketch). Because the shock is only slightly oblique, the 



shock relation concerning the product of velocity compo- 
nents normal to the shock front ahead of and behind the 
shock (ref. 7) may be written approximately as 

{V-U + ^,-u,)(V+^,-uO = 

2 \j (0.40.)] (la) 

The assumption of a slightly oblique shock also provides 
that the equations of conservation of normal energy and 
mass, respectively, may be written; 

~iV+^,-u,r+^{Q^+e,)= 

l(V-U+^-u,y+^.{e,+e,) (lb) 

J y — 1 

(R.-[-p.)(T>? -«.) = (7?2+P2)(T'-t/4« (Ic) 

The remaining oblique shock relation states that the velocity 
component parallel to the shock is unaltered by passage 
through the shock. Because the shock is assumed to be only 
slightly oblique, this relation is, approximately. 


Vl—V^y = V2— (F— f7)|y 

or, 

. V 2 — V 1 

(Id) 

The shock relations for the basic undisturbed shock propa- 
gation are obtained from equations (la), (lb), and (Ic), 
with small quantities neglected: 

(F C)F=2;^;(iF-+^F,e.) 

(2a) 


(2b) 

R,V=R2{V-U) 

(2c) 


Terms of equations (la), (lb), and (Ic) which are of first 
order in small quantities yield the conditions which the dis- 
turbance field must satisfy at the shock: 

( F- ^ [ F(| , - u,) + ^ 0,] 

(3a) 

V(tt — "a.)-!-— — ^0i = (F — U){^t — '*^ 2)4 -y02 (3b) 

y — 1 y — i 


— 'M'l) — P 2 (F — U) +7?2(?( — 142 ) (3c) 

Equations (2), the state equation 



Position of shock front if---'" '"--Instantaneous position of 

there were no disturbance disturbed shock front 


Figure 1. — Notation for shock wave propagating into a region of 
weak disturbance. 


and the assumption that the incident flow ahead of the shock 
is isentropic permit equations (3) to be simplified, yielding 
the following set of disturbance shock relations (eq. (Id) is 
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also included), which relate the disturbance fields ahead of 
and behind the shock: 


— W2 T> — "*^1 

F F 


(5a) 

II 

4-“ 

p Pi 

(5b) 

Pi 79 'U't t 

V ^ 

■ T) 7*' 

(5c) 

V2 vi U 
F^F F 

ky 

(5d) 


where the coefficients are constants depending on the Mach 
number of the undisturbed shock M—Vjai: 






7 + 

7—1 2 
7-(- 1 yA-I' 

9 




C, 


1 


7—1 


l+ll M-’ 


Dr- 


yAd 


2yAd ^ 

7—1 


Ad'^- 


2 

7 — 1 




2 7—1 


yAd 


( 6 ) 


EQUATIONS OF PLANE DISTURBANCE FIELD 

In addition to the perturbation shock relations (eqs. (5)) 
which are concerned with the compatibility at the shock of 
weak disturbance fields ahead of and behind the shock, the 
equations satisfied by the disturbance fields themselves are 
required. The nonviscous equations of motion are written 
in a coordinate system at rest relative to the gas ahead of or 
behind the shock. That is, the following equations apply 
in the coordinate system of figure 1 for the disturbance field 
ahead of the shock, and in a coordinate system moving with 
velocity U, for the disturbance behind the shock. Subject 
to interpretation of the coordinate system, the same equa- 
tions apply in both regions, and therefore the subscripts 1 
and 2 are omitted for the time being. To first order in small 
quantities: 


Momentum: 


Continuity: 

Energy: 

State: 




Vt- ^Vv 


P tArR{Ux-\-v-^ — 0 
c„Ret+P{Ux+Vy) = 0 

p_p f 

P~R'^Q 


( 7 ) 

(8) 

( 9 ) 

( 4 ) 


These equations may be combined to yield the equation for 
the pressure disturbance: 


, d^\n 2_7-P 


( 10 ) 


the equation for the entropy disturbance (which may be 
shown to be proportional to {pjP) — yipjR)): 





( 11 ) 


and the equation for the vorticity 
flow: 


bt 


{Vx—u„) = 0 


of the disturbance 

( 12 ) 


Thus, the pressui'e disturbance satisfies the wave equation 
(10), and any entropy variation (eq. (11)) or vorticity varia- 
tion (eq. (12)) is steady, relative to the main flow. Any 
disturbance field satisfying these linear equations may be 
regarded as composed of two parts, one steady and the other 
unsteady, in a coordinate system at rest in the main flow. 
From equations (11) and (12), variation of entropy and vor- 
ticity may be assigned to the steady flow; and from equation 
(8), the associated velocity components satisfy the incom- 
pressible continuity equation. From equations (7), pressure 
variations must be assigned to the unsteady flow, satisfying 
the wave equation (10). The unsteady portion of the flow 
may then be regarded as produced by a pattern of sound 
waves. A weak nonviscous disturbance field may therefore 
be considered to include: 

1. An unsteady, isen tropic, iiTotational disturbance, which 
ma\^ be regarded as produced by a pattern of sound waves, 
and 

2. A steady rotational disturbance of constant pressure 
and (in general) variable entropy and density. 

TYPES OF INITIAL DISTURBANCE CONSIDERED 

The present analysis concei'ns the interaction of a sliock 
wave with three types of initial plane disturbances: 

A. SOUND WAVE OVERTAKEN BY SHOCK 

The shock moves into a region in which a plane sound 
wave is propagating in a direction oblique to the dhection 
of shock propagation (fig. 2(a)). Since the shock velocity is 
supersonic relative to the gas ahead, it will overtake the 
sound wave, whatever its direction of propagation. The 
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solution for the interaction of such a wave with a shock may 
in principle be generalized by linear superposition to pro- 
vide analysis of the passage of a shock through any isen- 
tropic field of small disturbance. 

A general plane sound wave may be represented as follows 
(the particular profile of the wave need not be specified); 


where 


V 


= AJ 


V 

Pi 

Pl~ 

EL- 




V Xi J 


/ mxi—ly~\-ait'\ 


V Xi J 

•- 

l'mx\-ly-\-a\t\ 


V Xi J 


/ mxi — ly -\- ait~\ 


\ Xi J 

> 


Z=sin \ m=cos \pi 


(12a) 


(13) 


and Xi is a length characterizing the scale of the disturbance. 
If the function / were a sine wave, Xi would be equivalent 
to the wave length. By equations (12), (11), and (7), 
respectively, 

lAi= — mA2 


Pertinent equations of motion provide, as before, 

lAi = — mA2 '1 


Az=yA4, 


Ax= — m ^ As 


(14b) 


The coordinate system X 2 ,y is fixed relative to the flow behind 
the shock. 


C. STATIONARY VOBTICITY WAVE OVERTAKEN BY SHOCK 

The shock moves into a region occupied by a stationary 
plane shear disturbance of constant density, oblique relative 
to the shock front (fig. 2(c)). A system of such waves may 










(b) 


V-U 


xz 


As=yAi 


A: 


m 

yM 



(12b) 


The disturbance is of the type 1 discussed in the previous 
section (unsteady, isen tropic, irrotational), and is longi- 
tudinal ; that is, the fluctuating velocity component is in the 
dhection of propagation of the sound wave. 

B. SHOCK OVERTAKEN BY SOUND WAVE 

The sound wave propagates relative to the fluid behind 
the shock, in such a manner as to overtake the shock (fig. 
2(b)). Thus, consideration will be restricted to cases for 
which — (i 2 cos 1^21 — U. The initial disturbance may be 

specified in a manner similar to that employed for the pi-e- 
ceding case. 

The subscripts 1 and 2 have been introduced to denote the 
flow ahead of and behind the shock, respectively. In the 
present problem, the entire flmv disturbance occurs behind 
the shock. The subscript 2 is therefore appropriate to both 
the incident and reflected waves, which will hereinafter be 
distinguished by second subscripts 1 and 2, respectively. 


^ jY mx2—lyPa2t \ ' 
^ iV 1^21, / 

’^ = Asf + 

V V X 21 / 

Asf ( 

P 2 \ X21 / 


R 2 \ X 21 J J 


(14a) 



(a) Shock overtaking .sound wave. 

(b) Sound wave overtaking shock from behind, 
(cf Shock overtaking stationary shear wave. 

Fioure 2. — Type.s of initial disturbance considered. 


be employed to represent a turbulent field (refs. 3 and 4). 
Therefore, the effect of the passage of a shock through a 
single oblique shear wave may in principle bo generalized by 
Fourier superposition to provide an analysis of the passage 
of turbulence through a normal shock. 

The incident vorticity wave, of arbitraiy profile, may be 
i-ep resented as follows: 



From equation (9), 

mAi=lA2 (16) 

The disturbance is therefore a special case of the type 2 
discussed in the previous section (a steady vorticity disturb- 
ance of constant pressure), and, by continuity, must be 
transverse; that is, the fluctuating velocity component is 
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parallel to the plane of the shear waves. Since the wave is 
transverse, there may be a component of velocity disturbance 
parallel to the shock (perpendicular to the plane of fig. 2(b)) 
which may be of arbitrary amplitude. 

This type of interaction (problem C) is treated in refer- 
ence 6. 

ANALYSIS OF INTERACTION BETWEEN SHOCK AND 
IMPOSED DISTURBANCE 

THREE TYPES OP REFRACTION AT SHOCK 

It has previously been shown that the equations of motion 
imply that any weak disturbance field may be divided into 
two parts: an unsteady isentropic irrotational field, and a 
steady vorticity disturbance. This point of view may be 
adopted with regard to the disturbance downstream of the 
shock, produced by the interaction. 

Simple refracted sound wave — problem A. — In the case 
of a sound wave overtaken by a shock (problem A of the 
previous section), it would seem reasonable to expect that 
for near either 0 or tt (fig. 2(a)), the isentropic part of 
the downstream field would simply be a refracted sound 
wave traveling away from the shock. This is indicated by 
the sequence of events shown in figure 3(a). At time h, 
the initial wave intersects the shock front at point Pi. At a 
later time t 2 =ti-\-&t, the sound wave has moved a distance 
tti M, the shock has moved a distance V St, and the intemection 
occurs at point P 2 . In the meantime, a cylindrical sound 
wave has been generated at point Pi as a result of the shock 
interaction and expands with velocity U 2 , while being con- 
vected with a velocity U. Thus, at time ^ 2 , the effect of the 
intersection at h is felt within a cylinder of radius a 2 St, with 
center at point Qo. According to figure 3(a), an envelope is 
formed and may be identified as a simple refracted sound 
wave. 

Attenuating refracted pressure wave — problem A. — 

Figure 3(a) is drawn for a rather small value of 1 ^ 2 - If i'l is 
increased, there appears a critical value \pci (hg- 3(b)) be- 
yond which no envelope may be drawn. Thus, when ypi'^ipci, 
the influence of intersection Pi is felt at P 2 before the inter- 
section arrives at P 2 . However, as \f/i is further increased, 
there appears another critical angle \pcu beyond which simple 
envelopes may again be drawn (fig. 3(d)) and simple sound 
wave refraction occurs. 

Wlien the downstream pressure disturbance 

cannot be a simple sound wave. The cylindrical sound 
waves produced by the interaction at the shock do not 
coalesce, but rather continue to expand independently, thus 
diminishing in strength as time progresses. Accordingly, the 
isentropic part of the downstream disturbance may be ex- 
pected to die out at large distances downstream of the shock. 
This attenuating disturbance may, however, be expected to 
remain planar, because both the incoming disturbance and 
the shock are plane. This attenuating wave has been called 
a pressure wave rather than a sound wave, because, as will 
be shown subsequently, it does not propagate at the local 
velocity of sound. 

Expressions for Tpa and \j/cu hi terms of shock Mach number 


5 

may be obtained from the following equation derived by 
inspection of figure 3(b): 

02 ^- (F- C7)2= (F cot ^,-1-ai esc (17) 

The solution of this equation is shown in figure 4, labeled 
“sound wave.” The cimves labeled “stationary vorticity 
wave” will be discussed subsequently. The emves both 
approach a value of 180° for M—\, and have a half-order 
singularity there. As <» , the curves become symmetric 
about ^1 = 90°, because di becomes insignificant compared 
with F. The limiting value of rf/d is 67.8°. 

Steady vorticity wave — problem A. — If a vorticity disturb- 
ance is created at the shock-disturbance intersection (fig. 3) 
and is thence convected with a velocity U, a vorticity wave 
appears along the line cormecting P 2 and Q 2 , whatever the 
value of ^i. 

Thus, of the three types of refractions discussed, the 
“steady vorticity wave” always appears, in combination 
with either a “simple refracted sound wave,” or an “at- 
tenuating pressure wave”, depending on the angle ypi. 






(a) 

(b) = 

(c) cu- 

(d) lAcu 

Figure 3. — Formation of waves behind shock, because of interaction 
with sound wave (problem A). 


6 


EEPORT 1165 — ^NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


Expectations for problem B.— In the case of the shock 
overtaken by a sound wave, it may be shown in constructing 
figures similar to figure 3 that a simple reflected sound wave 
will always occur, in conjunction, of course, with a steady 
vorticity wave. 

Expectations for problem C. — In the case of the shock 
overtaking a stationary shear wave, sketches may be drawn 
similar to those presented in figure 3, except that the incident 
wave does not move in the time interval St. The qualitative 
character of the downstream disturbance is expected to be 
the same as for problem A. Despite the fact that the initial 
disturbance is not a sound wave, sound waves are produced 
by the interaction and form envelopes when or 

In this case, the values of and feu are ob- 
tained from equation (17) with the term involving Ui omitted. 
The solution, which is presented in figure 4, is symmetric 
about i/'i = 90°, approaches 90° at M=1 with a half-order 
singularity, and has the same asymptote at M= oo as in 
problem A. As in problem A, when an attenu- 

ating pressure wave occurs, and for all values of fi, the steady 
vorticity wave appears. 

The foregoing discussion ma 3 ^ bo summarized as follows: 
In case A, the incident sound wave refracts, as either a 
simple sound wave or as a more complicated attenuating 
pressure wave, and an oblique steady wave of variable 
entropy and vorticity appears. In case B, the incident 
sound wave reflects at the shock as a simple sound wave, 
and a steady vorticity wave appears. In case C, the initial 
“vorticity wave” refracts to form a stationary vorticity 
wave in which, because of the action of the shock, the entropy 
also varies. In addition, a sound wave, or pressure wave, is 
produced by the interaction. 

SOLUTION OF PROBLEM A 

In the problem of an oblique plane sound wave overtaken 
bj'' the shock, there are two different solutions to be ob- 
tained — one for and fcu<i'Pi<ir and another for 

The first is the simpler and will be presented 

first. 

Solution when 0C^fi<Cfci oi fcu<Cfi<^7r . — The initial dis- 
turbance ahead of the shock is described in equations (12) 
and (14). The discussion of the preceding section has es- 
tablished that, in these ranges of fi, the pressure variation 
behind the shock is associated with a simple sound wave. 
According!}^, the disturbance pressure is ivritten 

E|_A/(!Si+0±5L*) (18) 

where K, a, /3, and X 2 are to bo determined. In order that 
the pressure satisfy the wave equation (10), 

a" + )3^=l (19) 

Equation (18) is written on the assumption that the profile 
of the pressure disturbance carries through the shock un- 
distorted, though its orientation, magnitude, and scale may 
change. This assumption may be regarded as a trial, the 
correctness of which is inferred from the self-consistency of 
the entire solution so obtained. 


In view of the requirements of the shock relations (eqs. 
(5)), the arguments of the downstream pressure wave (eq. 
(18)) and the initial wave (eq. (12)) should match. At the 
shock, Xi= Vt and X 2 — (V~U)t. The matching requirement 
therefore is 


[a(F— ?7)-|-ffl2] t+l3y _ (rnV+a{)t—ly 
X2 Xi 

Equating coefficients of y and t yields 


( 20 ) 


where 


a 


1 

1 —r 




(8=-/ 


X2 

Xi 


I 


r=UIV 


( 21 ) 


( 22 ) 


Equations (19) and (21) yield a quadratic equation for 
X 2 /X 1 , the meaningful solution of which is 



where 



n=llm = ta,n 


(24) 



Fioube 4. — Critical angle of incident wave for formation of refracted 
sound wave behind shock. 
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The other solution corresponds to a wave moving in the 
same dhection as the shock and is rejected. (Inspection of 
fig. 3(a) shows that two families of envelopes might be indi- 
cated mathematically, and that only one is physically sig- 
nificant.) Thus, the inclination and scale of the pressm*e 
wave are fully determined, and only the magnitude K re- 
mains to he found. 

The vorticity wave will also be assumed to have a profile 
given by the function/. In view of equation (12), its argu- 
ment can be a function only of X 2 and y, and must further- 
more match the argument of the incident wave at the shock, 
in order that the shock relations may be satisfied. At the 
shock, the argument of the incident wave is given by the 
right-hand side of equation (20). Also, t=X 2 l(.V—U), at the 
shock. The argument of the vorticity wave must therefore 
be 


m-f 1/M 
1 — r 


X2—Iy 


Xi 


(25) 


The inclination and scale of the vorticity wave are thus de- 
termined. 

The density fluctuation behind the shock consists of two 
parts — one part associated with the pressure fluctuation to 
form the sound wave, its magnitude determined by eciuation 
(11); the other part associated with the vorticity wave. Thus, 


P2 y: / ax 2 -\-^y-\-g' 2 t \ ^ y 1 —r 

R2 y ' \ Xo / \ Xi 


(26) 


Likewise, the velocity components each consist of two 
parts, the first associated with the sound wave and the sec- 
ond associated with the vorticity wave. Accordingly, 




m- 


- 1 /M 


l—r 


-ly 


Xi 


(27) 


V \ ^2 / 




X2—Iy 


Xi 


(28) 


The requirement that all terms in the shock equations (5) 
have the same functional form and the same argument 
suggests that 


it=VLi r (^+l/-^ Vt-ly 
L Xi _ 




n 


‘l-fl/wiM 


^ j. r(m-l-l/M) Vt—ly 

X^ "J. 


(29) 


(Cross-differentiation shows that these two equations are 
compatible.) 

The solution is completed by the algebraic determination 
of the various unknotvn constants. The coefficients F and 
H may be found in terms of K through equations (7) ; / may 
be found in terms of G from the incompressible continuity 
equation. The remaining unknowns K, Q, O, and L may be 
successively determined by use of the four shock relations 


29H189—H tZ. 


(eqs. (5)) when it is recalled that the arguments of aU quan- 
tities have been arranged to match at the shock. Neither 
the details of the remaining procedure nor of the final solu- 
tion are particularly interesting, and therefore the analysis 
has been completed in appendix B. The numerical results 
will be discussed in a subsequent section. 

Form of attenuating pressure wave when/'d<C'Ai<('Ai!«' — It 
is intended to form the solution for this range of in essen- 
tially the same way as was done in the preceding paragraph 
An essential step in that solution was the assumption tha 
the refracted sound wave has the same profile as the initial 
disturbance. Therefore, in order to proceed with the analysis 
of the case it is first necessary to determine the 

form of the pressure disturbance just behind the shock and 
the manner in which it attenuates with distance behind the 
shock. 

Tentatively, the pressure is written as a function of two 
variables only; 

in, f) (30) 

where 

, r 

f iax2 + ^y + cVt) 

Al J 

and the constants d, a, j3, and c require determination. The 
variable 17 is proportional to distance behind the shock front, 
—X 2 -\-iV~U)t. Neither the undisturbed shock front nor 
the incident wave has curvature. Therefore, it is expected 
that along any one of a family of planes moving with con- 
stant velocity, any variation of pressure would be due solely 
to the attenuation associated with distance behind the shock. 
This consideration leads to the definition of the second vari- 
able f, such that the equation constant defines a plane 
moving obliquely with a constant velocity. 


The wave equation (10) is satisfied if 



(32 a) 

a=——2 (l-^)c 

d2 

(32b) 


Z! 

0.2^ 


[d^(l -ry-F] = d^-ia^+ 13^ 


(32c) 


The boundary conditions to be applied in solving La- 
place’s equation (32a) are 


#(“,f) = 0 (33) 

^2 

and a condition (eq. (5c)) at the shock providing compati- 
bility with the initial disturbance. In the following discus- 
sion, this information will be used to infer a Iffiely form for 
the pressure wave. 

Part of the downstream velocity variation is associated 
with the pressure to form an isentropic irrotational flow. 
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For this part of the disturbance field, a velocity potential 
<p(ri,^) may therefore be defined such that 


d . a 

V2 = <Py = Z- <Pt 


(34a) 

(34b) 


In view of equations (7) and (34), 


1 


^ grad <P=-^ gi'ad po_; 


dt 


. Pk 


Rt 


-<Pt 


Therefore, 




P2 


(34c) 


In the case of the simple sound wave, compatibility at the 
shock was obtained by supposing that the profile of the dis- 
turbance (the function/) carried through the shock undis- 
torted. In the present case, the corresponding assumption 
would be that, just behind the shock at t; = 0, the various 
disturbance quantities each contain a term proportional to 
jit) and that t should therefore match the argument of the 
initial disturbance at the shock. From equation (34b), if 
^ 2 ( 0 , t) is to contain a term proportional to j(t), then the 
potential must contain a term proportional to 

= r fir) tan-i dr 
7T J-= 7? 


which satisfies equation (32a) and hence the wave equation, 
satisfies boundary condition (33), and has the property that 
(0, (■)=/(?). This solution may be regarded as the result 
of a distribution of singularities along the plane of the shock 
( 77 = 0 ). In the skew coordinate system 77 , C, these singulari- 
ties may be identified as potential-flow vortices. Therefore, 
from equations (34a) and (34c), and p 2 would each contain 
a term linear in fit) at the shock and, in addition, a term 
linear in 

5(r)^^/«(0,f)=ip.F. (35a) 

7T T— f 


where P.V. denotes the Cauchy principal value of an im- 
proper integral. 

Of course, V 2 would lilcely contain a term linear in ^(f) at 
the shock also, and therefore the potential would have an- 
other part 

= r fir)ln[r,^+ir-mdr 

Zi'K J — CO 


wnere 




dr 


/ ( 37 ) 


✓ 

and the constants KP-'' and require determination. At 
the shock, equations (35) and (36) provide that: 

2^^=if»)/(f)+if(^>5(f) 

Examples. — (a) If it happens that /(f) = sin 27rf, then 
,?(f)=cos 27 t^, and 


Thus, 


$(“ = 6 sin 2-wt 

= cos 2irt 

( 1 ) gjj^ cos 2xf] 

Jr 2 


(38) 


and, therefore, the disturbance undergoes a phase shift in 
passing thi’ough the shock and subsequently attenuates ex- 
ponentially with distance behind the shock. 

(b) If /(r) = (l + H~‘, then ^(r) = — + (see ac- 
companying sketch), and 



Solution when ypci^4'i<^'Pcu - — The form of the pressure 
wave has been adduced in the previous paragraph (eq. (36)). 
The quantities d, a, fi, and c may be found by using equations 
(32) and the requirement that the disturbance function / 
have the same argument ahead of and behind the shock. 
At the shock, Xi=Vt, X 2 — (C— U)t, this requirement leads to 
the equations 

m-f 1/M— c 

a = 

1— ?■ 

Equations (32b) and (39b) yield 


(39a) 

(39b) 


satisfying equations (32 a) and (33), and having the proper- 
ties: 

w'” (0,f) =git ) ; (0,f) == -fit) (35b) 

This solution may be regarded as a distribution of potential- 
flow sources along the plane of the shock in the n,t system. 

Thus, the quantities associated ivith the attenuating pres- 
sme wave may tentatively be written in the following form: 


P2iv,t)_ 


m-|- 1/M 


.Z! 


(1-r)^ 


and equation (32c) may be solved for d: 


^ 






(39c) 


(39d) 


:/<-<‘>$™(77,f)-bif«$'2)(^,f) 


(36) 


The vorticity wave is expected to involve a linear combi- 
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nation of the profile functions / and g, just behind the shock, 
and is not expected to change its form subsequently, because 
it is not time-dependent in the X 2 ,i/ coordinate system and 
therefore cannot attentuate, as does the pressure wave, or 
otherwise change character. A matching procedure at the 
shock yields, as previously, the argument given in expres- 
sion (25). 

The remaining analysis parallels that following expression 
(25) in the paragraph Solution when or 


{[aV(l-r) + a 2 ]t+^y]^~ ilmV(l-r) + a 2 jt-ly} (45) 
A 22 ^21 

whence, 

(46b) 

A 2 I 


rv2 7-^2 




/m+lIM , \ 




+ 






(40) 


(41) 


Equations (19) and (46) provide a quadratic equation for 
1 ^ 22 A 21 , the useful solution being: 


X 22 

X 21 


(I 2 


-H-ry 


( 1 - 


-r)"d-2 ||m(l— r)-t-|2_ 


(46c) 


The other solution is X 22 /X 2 i=l, corresponding to the incident 
wave itself. 

The right side of equation (45) is used as the argument in 
the expressions for and 




^m-f 1/M 


Xi—ly 


+ 


X. 


7 (2,^ 


( 


1—r 


X2—Iy^ 


|~ (m+l/yFt-Zy J ^ ^ 


{m -\-llM)Vt—ly 


Xi 


(42) 
■] (43) 




1 1 




m.-\- llM)Vt — ly 

In 




+ 


^,= VLf{[mV{l-r) + a2]t-ly} 




ffl2 

mV 


Lf{[mV{\—r) + a2]t—ly} 


H47) 


The right side of equation (45) also yields the argument of 
the vorticity wave, when the substitution t=X 2 j{V — U) is 
made: 


X 21 


(48) 


L (2) g |^ (m + l/M)Ff- J y j| 

The various unknown constants remain to be determined 
algebraically through equations (7), the incompressible con- 
tinuity equation, and the shock equations (5), as before. 
These equations suffice to determine a greater number of 
constants than were required in the previous case because 
the functions/, g, and and their derivatives form 
two separate groups of functions whose coefficients may be 
separately equated. Details of this procedure are provided 
in appendix B. 

SOLUTION OF PROBLEM B 

The analysis of problem B (sound wave overtaking shock 
from behind) is identical in all essential respects to that of 
problem A when only a simple refracted sound wave is in- 
volved. The only differences which arise are the slightly 
different matching of arguments at the shock and a slightly 
different form assumed by the shock relations (5) . 

The equation (18) is adopted in the present case to rep- 
resent the reflected sound wave. Matching the argument of 
this expression with that of the initial disturbance (eq. (14)) 
yields 


Equations (26), (27), and (28) may be adopted to complete 
the description of the flow, except that expression (48) must 
be used for the argument of the vorticity wave. The analysis 
is completed in appendix B. 

SOLUTION OF PROBLEM C 

The only differences between problem A and problem C 
(stationaiy vorticity wave overtaken by shock) involve the 
shock equations (5) and the matching of arguments at the 
shock. The difference in matching is due to the fact that in 
the present case, the disturbance is stationary, while in prob- 
lem A, the disturbance moves with velocity a,. Accordingly, 
when 0<C'f'i<ifci, equations (18), (21), (23), and (25) through 
(29) apply directly to the present case if the quantity 1/M 
is omitted wherever it appears explicitly. When 
equations (36), (37), and (39) through (44) may also be 
adopted, again by omitting terms proportional to 1/M. The 
remaining details of the analysis are provided in appendix B. 

As previously mentioned, an initial disturbance of type C 
may have a third fluctuating velocity component parallel to 
the shock, which might be represented as follows: 
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(see eqs. (15)). The amplitude Ag is arbitrary, within the 
limitations of linear analysis. This disturbance will pass 
through the shock unaffected and become part of the steady 
vorticity wave behind the shock. Thus, 

The shock front itself will not be affected by this purely 
transverse disturbance. 

SINGULARITIES AT AND 

Equation (23), which applies in the refracted sound-wave 
solution, contains a radical which vanishes at if'ci and 
and becomes imaginary when tpciK'PiKi'cu- It may be shown 
that the quantity under the radical vanishes with a nonzero 
slope. Accordingly, the quantity X 2 Ai, though finite, has a 
half-order singularity at ^l/ci and \pcu- Furthermore, this 
quantitj^ is involved in all the formulas characterizing the 
refracted sound wave (see appendix B) . In the range \}/ci < 
'Pi<'Pcu, the attenuation coefficient d (eq. (39d)) vanishes 
with half-order singularities at ipci and tpcu and similarly 
affects the remainder of the analysis. 

The reason for this singular behavior ma}’ be inferred from 
figure 3. According to figure 3(a), when the 

refracted sound wave is the envelope of an imagined succes- 
sion of cylindrical waves, as shown in the following sketch. 



At i/'c!, however, according to figure 3(b), the cylindrical 
waves all meet at a comm.on point of tangency, as shown in 
the sketch. Therefore, successive waves reinforce at one 
point, giving a singularitj^ of the flow. This singularity, of 
course, depends on the fact that the theory is linear. An 
exact anal 3 ^sis would presumably show steep, though not 
singular, flow gradients. 

This situation is similar to that arising in the linearized 
analysis of compressible flow about bodies: as the i\Iach 
number of the flow approaches 1 , the i\Iach waves have a com- 
mon point of tangenc 3 " at the nose, and the wave drag shows 
a reciprocal half-order singularit.y in Mach number. In the 
present case, the physical quantities remain finite, but have 
infinite rates of variation with 4/i. 

The preceding discussion applies to problems A and C, 
but not to B. 

RESULTS AND DISCUSSION 

In the following paragraphs, the results (presented in 
graphical form) will be described for problems A and B. 


Results for problem C are presented in reference 6 . The solu- 
tion for each problem has essentially three elements: 

1 . Disturbance of the shape of the shock front 

2 . Characteristics of the isentropic pressure wave behind 
the shock 

3. Characteristics of the steady vorticity wave behind the 
shock. 

Computations have been carried out for three Mach num- 
bers, 1 , 1.5, and infinity. Of course, the case M=1 is really 
degenerate, because the shock is then a weak sound wave, 
and therefore the interaction with the incident disturbance 
is obtained by linear superposition, the initial disturbance 
passing through the “shock” with no change. Also, Af=l 
is a singular point, because the range of angles of incidence 
of initial disturbance for which the attenuating pressure 
wave appears in problems A and C vanishes with a half -order 
singularity (fig. 4), and in problem B, because the incident 
wave is unable to overtake the “shock” moving with sonic 
velocity. 

The results will show’’ that the critical angles and 
are also singular points. In many instances there are not a 
sufficient nu.mber of points near the singular points for which 
computations have been made, so that not aU the curves can 
be faired with complete confidence. For this reason, the 
computed points are showoi circled, so that the basis for the 
fairing will bo clear in each case. 

No mention will be made of the temperature disturbance 
behind the shock, which maj’’ be obtained directly from equa- 
tion (4) if the pressure and density disturbances are known. 

All disturbance quantities found by the linearized analysis 
of the present report will be proportional to the intensity of 
the incident wave. Therefore, results are divided by the 
pressure amplitude of the incident w^ave A 3 (see eqs, ( 12 a) 
and (14a)). 

PROBLEM A— SHOCK OVERTAKING SOUND WAVE 

1. Shock front disturbance. — In figure 5 are shown the 
amplitudes and ' =i®/A 3 of the incremental 

velocity of the shock front, due to the interaction. From 
equations (29), (43), and (44) these amplitudes are associ- 
ated, respectively', with the functions / (w^hich defines the 
profile shape of the incident wave) and g (which is an addi- 
tional profile function arising wdien 4'ci<^'Pi<CPat) to give the 
actual incremental velocity. The variations with xpi are 
quite extreme, particularly at the critical angles, where, in 
fact, there are half-order singularities. The variations with 
Mach number are equally^ severe. 

When \pi = 0 (incident wave moving parallel to and toward 
the shock), the figure shows that the shock front is retarded 
by a pressure wave. In the case M=l, this is because the 
velocity in an incident compressive sound wave, relative to 
which the shock (really^ a weak compression wave) propa- 
gates as it moves through the disturbance, is directed against 
the shock. When the incident compression wave moves in 
the same direction as the shock {\pi = ir), the shock front is 
speeded up for low shock IMach numbers; at i\/=l, this is 
true because the incremental velocity due to the incident 
wave is in the same direction as the “shock” movement. 
Whether \pi = 0 or w, there is a smaller accelerating effect due 
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Figure 5 . — Problem A; shock front disturbance (shock overtaking 

sound wave); 


to the higher velocity of sound in the incident compression 
wave — that is why the curve for M=1 is not perfectly 
antisymmetrical about i/'i=x/2. 

When M= , the curves are symmetrical about i^i = ir/2 
because the incremental sound-wave velocity is vanishingly 
small compared with the shock velocity. 

For each value of M there is a value of yp\ for which the 
shock intersects the incident sound wave permanently at 
one point on the traveling wave, and the problem becomes 
essentially steady, so that the increment in shock velocity 
vanishes (though a steady displacement occurs). This hap- 
pens when ai=—mV or ^i=sec“'(— M), yielding 131.8° 
when M=1.5, and 90° when M=°°, and is the case of 



Figure 6. — Problem A: inclination of refracted pre.ssure front (shock 
overtaking sound wave). 


steadji interaction of a Mach wave and a normal shock which 
has been treated bj^ Adams (ref. 2). 

2. Characteristics of pressure wave. — (a) Inclination of 
refracted pressure front: Figure 6 shows the angle between 
the directions of propagation of the shock and the refracted 
pressure wave behind the shock. In view of equations (18) 
and (31), this quantity is given by the equation 

i/'2 = cot"*( — “/(3) 

Of course, when this wave is not a sound wave, 

and the inclination shown refers to a front parallel to which 
physical quantities depend only on distance behind the shock. 
Outside the range the pressure wave is a sound 

wave. When M= , the curve is symmetrical about ^i = 7t/2. 





imiw 





60 70 80 90 100 110 120 130 140 150 160 

Inclination of incident wove, i//|, deg 

Figuee 7. — Problem A: attenuation of pressure wave behind shock 
(shock overtaking sound wave) . 

(b) CoefHcient of attenuation of pressure wave: Figure 7 
shows the coefficient d which appears in equations (31) and 
(39d). This quantity vanishes at \pci and indicating that 
the sound-wave solution and attenuating wave solution meet 
continuously at the critical angles. When M is either 1.5 or 
CO, the maximum value of d is about 1. This implies a 
rather rapid attenuation — if /(f) = sin 2xf, it has been shown 
(eqs. (38)) that the wave attenuates as exp{—2Tni) . From 
the definition of j? (eq. (31)), when d«l, the attenuation 
factor becomes at a distance behind the shock approxi- 
mately equal to l/2x times the wave length of the incident 
sound wave. 

(c) Propagation velocity of pressure wave when 

\pcu- Figure 8 shows the quantity c of equations (31) and 
(39c), combined with other quantities to give propagation 
velocity as a fraction of the speed of sound a^. Since the 
solution when meets the sound-wave solution 

continuously, the propagation velocity is equal to the speed 
of sound at ^pei and \pcu- The change of sign of c is taken 
into account in figure 6 by the 180° shift of direction shown 
at the angle for which c = 0. 

(d) Katio of scales of pressure waves behind and ahead of 
shock: Figure 9 shows the quantity X 2 /Xi (eq. (23)), which 
was defined only for the sound-wave solution. However, 
inspection of equation (31) shows that the equivalent quan- 
tity when ^ci<'Pi<'Pcu is 

The reversal of sign of X 2 /X 1 signifies a reversal of the direc- 
tion of propagation of the pressure wave relative to the shape 
of the incident wave. For example, when i^i = 0, the incoming 
and outgoing waves might appear as follows: 

I ^ - Incident-wave profile 


-X2A1 positive 



100 no 120 130 140 150 160 

Inclinotion of incident wove, deg 


170 180 


Figure 8. — Problem A: propagation velocity of pressure wave behind 
shock (shock overtaking sound wave). 
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Figure 9. — Problem A; ratio of scales of pressure wave upstream and 
downstream of shock (shook overtaking sound wave). 


whereas when iAi = ir, they would appear thus: 


Xj/Xj negofive 


The difference between these two cases consists of a difference 
in sign in the arguments of the refracted wave in the two 
cases, and arises formally in the present analysis as a change 
in sign of X 2 /X 1 . 


Shock 
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Inclination of incident wave, »//j, deg 

Figure 10. — Problem A: amplitudes of pressure disturbance behind 
shock (shock overtaking sound wave). 


When M—1.5, the magnitude of X 2 /X 1 is greater for i/'i = ir 
than for t^'i== 0 , because, when the incident wave is traveling 
in the same direction as the shock, the shock requires a 
longer time to traverse the incident wave than when the two 
waves travel in opposing directions. 

(e) Amplitudes of pressure disturbance behind shock; 
Figure 10 shows the coefficients (of /, or of and 
when which describe the pressure wave behind 

the shock: (see eqs. (18) and (36)). As in 

figure 5 , the flow is shown to be singular at the two critical 
angles and to vary markedly with both ^1 and M. At 
M=1.5, near iAi = e, the refracted sound wave is seen to be 
very weak. 



Figure 11. — Problem A: coefficient for part of uz associated with 
pressure wave (shock overtaking sound wave) . 


Since the pressure wave is isentropic, the coefficients of 
the corresponding part of the density variation are equal to 
1/7 times the pressure coefficient. 

(f) and (g) Coefficients for velocity components in pres- 
sure wave: Figures 11 and 12 show the coefficients for the 
velocity components associated with the isentropic pressure 
wave behind the shock (eqs. (27), (28), (41), and (42)). The 
velocity resultant is longitudinal with respect to the direc- 
tion of propagation of the pressure wave, except when 

'f'ciK'f'iK'P CU* 
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Figure 12. — Problem A: coefficient for part of Vo associated with pres- 
sure wave (shock overtaking sound wave). 


3. Characteristics of steady vorticity wave behind shock. — 
(a) Inclination of the steady vorticity wave; Figure 13 
shows the inclination of the vorticity front behind the 
shock. From equation (25), 

l/-3 = COt“' 

(b) Ratio of scale of vorticity wave to that of incident 
sound wave: From equations (12a) and (25), this is 


-\-llmM~ 

n{l—r) 



Figure 13. — Problem A: inclination of vorticity wave behind shock 
(shock overtaking sound wave). 


5 c 


O C 

or o 



Inclination of incident wove, «/^ 2 > 


Figure 14. — Problem A: ratio of scales of vorticity wave behind 
shock to that of incident sound wave (shock overtaking sound 
wave) . 


which is shown in figure 14. When Ad— I and = the 
scale ratio goes to «> because the “shock” is unable to over- 
take the incident wave. 
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PifiUKE 15. — Problem A: amplitude of dcn.sity variation in vortieity 
wave behind shock (shock overtaking sound wave). 
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Figure 16. — Problem A: coefficient for part of associated with vor- 
ticity_wave (shock overtaking sound wave). 
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Figure 17. — Problem A: coefficients for part of vn associated with 
vortieity wave (shock overtaking sound wave). 



Figure 18. — Problem B: shock-front disturbance (sound wave over- 
taking shock) . 


(c) Amplitude of density variation in vortieity wave be- 
hind shock: Figure 15 shows and Q®' (eqs. (26) and 

(40)). Apparently, when M is of order 1.5, the vortieity 
wave is very weak. 

(d) and (e) Coefficients for velocity components in vor- 
ticity wave: Figures 16 and 17 show the coefficients for the 
transverse velocity field associated with the vortieity wave. 
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Figure 19. — Problem B: inclination of reflected sound wave (sound 
wave overtaking shock). 


PROBLEM B— SOUND WAVE OVERTAKING SHOCK FROM BEHIND 

Consistent with the previous analysis, computations have 
been carried out only for values of i> 2 i (incident-wave incli- 
nation) sufficiently close to 180° that the component of 
propagation velocity in the direction of the motion of the 
shock is greater than the velocity of the shock relative to 
the fluid behind. These values of are 134.5° and 112.2° 
for M=1.5 and °=, respectively. 

1. Shock-front disturbance. — Figure 18 shows the ampli- 
tude of the incremental shock-front velocity L'=LIA^. 
When ^21 is near 180°, an incident pressure wave displaces 
the shock ahead. 

2. Characteristics of reflected sound wave. — -The down- 
stream pressure wave in problem B is always a simple sound 
wave. 

(a) Inclination of reflected wave; The inclination 
cot“' (— a//3) (see eq. (18)) is shown in figure 19. In effect, 
the incident and reflected waves coalesce into a single wave 
at the critical angles. 

(b) Ratio of scales of reflected and incident sound waves; 
This ratio is given by equation (46c) and is shown in figure 
20. At i/' 2 i = x, the scale ratio is greater when M=1.5 than 
when M= co , just as in problem A, and for the same reason. 
At the critical angles, the ratio becomes 1 , because the incident 
and reflected waves coalesce. 

(c) Amplitude of reflected sound wave; The pressure 



Figure 20. — Problem B; ratio of scales of reflected and incident sound 
waves (sound wave overtaking shock) . 



Inclination of incident wave, ^g| , deg 


Figure 21. — Problem B: amplitude of reflected sound wave (sound 
wave overtaking shock) . 

amplitude K'^KjA^ (eq. (18)) is shown in figure 21. At 
i/' 2 i = ir, an incident compression wave reflects as an expansion 
wave with a strength which is greater the higher the Mach 
number, but always less than that of the incident wave. 
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Figure 22. — Problem B: inclination of vorticity wave produced 
behind shock (sound wave overtaking shock). 


The velocity components associated with the reflected 
sound wave are obtained simply from K' and are therefore 
not plotted. 

3. Characteristics of steady vorticity wave. — The coeffi- 
cient Q' of the density fluctuation in the vorticity wave is 
simply proportional to L' (combining eqs. (B26) and (B27) 
of appendix B) and is therefore not shown in a figure. 

At M=1.5, Q'=-0.147i' 

At M= cx=, Q'=-1.43i' 

(a) Inclination of vorticity wave: From equation (47), the 
inclination is given by 

^,3=cot-|^[l + ^^|f^]| 
which is plotted in figure 22. 

(b) Ratio of scale of vorticity wave to that of incident 
sound wave: As in problem A, this ratio is given by 
(1/Z) sin i /'23 and is shown in figure 23. As for the reflected 
sound wave, the scale ratio is larger for M nearer one. 

(c) and (d) Velocity variations in vorticity wave: Coeffi- 
cients G' and /' of the transverse velocity fluctuation (eqs. 
(27) and (28)) in the vorticity wave are shown in figures 24 
and 25. 

4. Wave reflection at critical angle. — The analysis and 
figures show that at the critical angle the incident and re- 
flected sound waves coalesce to form a single sound wave. 
This statement may be interpreted to mean that a sound 
wave incident at the critical angle reflects as a steady vor- 
ticity wave only. Therefore, the shock disturbance and 
vorticity wave characteristics may be expressed in terms of 
the pressure amplitude of a single incident sound wave of 
strength As-\rK=As{l-\-K'). 
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Figure 23. — Problem B: ratio of scales of vorticity wave behind shock 
and incident sound wave (sound wave overtaking shock) . 



Figure 24. — Problem B: coefficient for part of M 22 associated with 
vorticity wave (sound wave overtaking shook). 
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Figure 25. — Problem B: coefficient for part of Vzz associated with 
vorticity wave (.sound wave overtaking shock). 
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CONCLUDING REMARKS 

In principle, the interaction of a shock with any weak 
flow fleld may be obtained by first constructing the initial 
flow field as a linear combination of plane waves of varying 
strength and orientation. From the present analysis the 
interaction of each constituent wave ivith the shock may be 
found. Assembling the resulting waves behind the shock 
then would yield the desired solution. Regrettably, the 
formulas for the interaction depend on the angle of incidence 
in a rather complicated way and it would in general be diffi- 
cult to evaluate explicitly integrals in which these formulas 
are used for the distribution functions. Numerical pro- 
cedures could be used for this purpose, though a technique 
would be required for dealing with the singularities at the 
critical angles \pci and 

The nature of the solution of problem B perhaps requires 
a clarification, in that the angles of incidence yj/^i havejbeen 


restricted to the range for which the soimd wave overtakes 
the shook. For purposes of superposition, however, all inci- 
dence angles must be considered. Inspection of figures 6 to 
12 shows that a sound wave of any incidence angles between 
0 and 7T may bo identified either as an incident wave or as a 
reflected wave, in the sense of the analysis. For purposes 
of linear superposition, the distinction between incident and 
reflected waves is of no significance. The point of view may 
be adopted that whatever its angle, a constituent wave of 
the initial flow has associated with it another sound wave of 
another angle and a steady vorticity wave. The notion of 
cause and effect is not needed. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, February 8, 1964. 


APPENDIX A 

NOTATION 


The following symbols are used in this report: 

A], Ao, Aj, A,i coefficients of incident disturbance (eqs. 

(12a), (14a), or (15)) 
a velocity of sound 

coefficients in shock relations (eqs. (6)) 
coefficients in shock relations (eqs. (6)) 
dimensionless velocity of propagation of 
pressure wave behind shock (eq. (31) or 
(39c)) 

specific heat at constant volume 
coefficients in shock relations (eqs. (6)) 
coefficient of attenuation of pressure wave 
behind shock (eq. (31) or (39d)) 
coefficient of part of associated with pres- 
sure wave (eq. (27) or (41)) 
profile function of incident wave (eqs. (12a), 
(14a), or (15)) 

coefficient of part of associated with vor- 
ticity wave (eq. (27) or (41)) 
additional profile function appearing behind 
shock 


Bu B, 

Cl, 62 


c, 

Di, D2 

d 

F 


J 


G 


g 






-, eq. (35a)^ 


Q 

R 

r 

t 

U 

u 

T’ 


V 


Xu X-2 


1 / 

a, 

T 


H 


I 

J 

K 

L 

I 

M 

m 

n 

P 

P 


coefficient of part of V 2 associated with pres- 
sure wave (eq. (28) or (42)) 
functions involved in solution when <f'ci.<i 
i'lK'Pcu (eqs. (B19), appendix B) 
coefficient of part of V 2 associated with vor- 
ticity wave (eq. (28) or (42)) 
gas constant (eq. (4)) 
coefficient of p 2 (eq. (18) or (36)) 
coefficient for ^(y,t) (eqs. (29) or (43)) 
sin xpi or sin 4^21 

Mach number of shock (=T7«i) 
cos or cos i/'2i 
tan t/'i or tan ip 2 i 
mean static pressure 
perturbation in static pressure 


ri 

0 

e 

X 

p 

<x, X 
<!> 

<P 


4'h 4'2\ 


coefficient for density fluctuation in vor- 
ticity wave (eq. (26) or (40)) 
mean gas density 
ratio of U to V 
time 

mean velocity of gas behind shock (fig. 1) 
perturbation of velocity component in x- 
direction (fig. 1) 

mean velocity of propagation of shock in 
gas at rest (fig. 1) 

perturbation of velocity component in y- 
direction 

coordinates measured in the direction of the 
shock propagation, relative to which the 
gas is (on the average) at rest ahead of 
and behind the shock, respectively (figs. 
1 and 2) 

coordinate orthogonal to Xi or X 2 (figs. 1 and 

2 ) ^ 

functions defining pressure front (eqs. (21) 
or (39)) 

ratio of specific heats ( « 1 .4 for air) 
variable upon which pressure wave depends 
when 4'ci<4'i<4'cu (eqs. (31)) 
variable upon which pressure wave depends 
when 4'oi<4'i<4'cu (eqs. (31)) 

.mean static temperature 

fluctuation in static temperature 

scale of pressure wave (eq. (23)) 

displacement of shock front (fig. 1) 

fluctuation in gas densit3^ 

functions appearing in solution for L (eqs. 

(B9) or (B3I), appendix B) 
function associated with pressure wave when 
fciK'f'iK'f'cu (eqs. (37)) 
velocit}^ potential associated vuth pressure 
wave when 

angles of inclination of incident waves (fig. 

2) 
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^ 2 , 'p 22 angles of inclination of pressiure wave be- 

hind shock (figs. 6 or 19) 

^ 3 , 4'23 angles of inclination of vorticity wave behind 

shock (figs. 13 or 22) 

'Pcu lower and upper bounds, respectively, of 

the range of for which the attenuating 
pressui'e wave appears (eq. (17)) 

- Subscripts ; 

Subscript notation for partial differentiation has been 
used where convenient. 

1 conditions ahead of shock 


2 conditions behind shock 

Double subscripts 21 and 22 (problem B), incident and 
reflected waves, respectively. 

Superscripts : 

(1) coefficients associated with/ at shock 

(2) coefficients associated with g at shock 

Primed coefficients are referre d to inten sity of incident 
wave. (As in problems A and in problem C) . 

Primes are also used to denote ordinary differentiation of 
/ and g with respect to their arguments. 


PROBLEM A 

The solution of problem A is found in two separate ranges 
of ^ 1 , for which the solution involves a refracted sound wave 
in one case and a refracted attenuating pressure wave in the 
other. The analysis will be completed in that order. 

1. — Equations (27) and (28) are 

substituted into equations (7), and coefficients of f are 
equated, yielding: 


APPENDIX B 

COMPLETION OF INTERACTION ANALYSIS 

The shock equations remain to be satisfied ; it should be noted 


F=—%Koc 
7 V 


(Bl) 


When equations (26), (27), and (28) are substituted into 
continuity equation (8), the quantities associated with the 
sound wave combine to satisfy continuity. The terms asso- 
ciated with the vorticity wave then must satisfy the incom- 
pressible continuity equation (Ux-\-Vy= 0 ) because the corre- 
sponding density term is time-independent. Therefore, 
equating coefficients of/' yields 


_1 + 1/Mm 


n (1 —r) 


G 


(B2) 

1 


that the arguments of /for the various quantities have been 
matched at the shock, so that coefficients of/may be equated. 
From equations (5a), (27), (29), and (12a), 


G=L-F~By(L-A,)- B.,As 


(B3) 


From equations (5b), (5c), (5d), respectively, with the neces- 
sary substitutions made, 

Q = C\(L — Add-ftAs — — (B4) 


K — (L — Ai) -\- D 2 A 3 

^_ 1-|-1/Mm 


nr 


(H+I-A 2 ) 


(B5) 

(B6) 


'I'he terms H and / in equation (B6) maj^ be expressed in 
terms of L by using equations (Bl), (B2), (B3), and (B5). 
The resulting equation for L jdelds the result 


L'^ 


L 


CL2 

' 7 v 




yM 


A 3 


(B7) 


rn 


0/2 


1 -h 


yV 




n(l 


Mm 


Mm/ as n 


yV 


i-Pi) 


or, from equations (6), (21), and (23): 


m 






1 M /3 — Y 1 y 


4 m \7-(-l 




where 


and 


yM 


n(\ — r) 




(B8) 


1 + 


( ^ 


1 

Mm 

l + <r^ 
-1 


1- 


) l/2 


(B9) 


and where use is made of the results from equations (2), (6), 


and (22) that 


y+ 


1 (* w) 




When L' is known, the remaining coefficients follow from 
equations (Bl) through (B5). 
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2. — Equations (36) and (41) substituted into 

(7) give 
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^»>=A(i"-^i)+A^3 


«) a)$^a) (2>$^^(2>) (BIO) 

From equations (31) and (37), 


o. « = ( F- C7) ;r 

Ai Ai 

Al Ai 


( 1 , 2 ) 


Equations (37) shoAV that 




(2) 


Accordingly, equation (BlO) may be written; 

(1 -r)dF'“$/» +(1 -r)dF®$/^' + cF®$,«> 


= ~%5 [- + o^C®4>,<‘>] 

■y K 

By equating coefficients of and separately. 


(1 -r)dF<» +cF® =-%-JdK^» - oK^^) 
y V 


V (BID 


( 1 - r) dF - cF (di?C »> + ) 

Similarly, equations (36), (42), and (7) yield: 


(l-DdH^ + cA'^ 


(T2 


(l-r)dH®-cH<»=-^/3F:'” 

7 K 


7^2 

2 


/3if 


(2) 


(B12) 


When equations (41) and (42) are substituted into equation 
(8), and coefficients of/' and g' are equated, the following 
equations are obtained: 


n(l — r) 

1+^r^ 

J (2) = Q (2) 

n(l — r) 


(B13) 


} 


XO) = ^ /jjr(i) ^ J(l) ) 

nr ' 

j^(2) ^ (jji2) -|-7®) 


nr 


(B16) 


(B17) 


Equations (Bll) through (Bl7) may be combined to yield 
the solutions for and Z® 


Z<»' = ^ - 


hj\Jh 


( rn Z> 2 \ 

' VtM+A/' 


A, 






(BIS) 


where 


A.— 

/l3= 


m 

>' O’ 

l-r 1-A 
r v' 

2/(1 -r) 


A+A - 






o+vo 


(B19) 


Equations (B16) are substituted into equations (Bll), 
which may be solved to yield 


F®' = /iA®' + A(Z“>' I I -^2 


yM ' D 


) 


(B20) 


where 


hd — 




T+1 (l + <r^) 
2rdc{l-r) 


(B21) 


The shock conditions remain to be applied. Equations 
(5a), (5b), (5c), and (5d) yield, respectively, the following 
four pairs of equations, when the coefficients of / and g are 
separately equated: 


-A(Z® - AO-Z2A3 
(?® = Z®-F®-ZiZ® 

= A(Z <« - AO + A Aa- A«>/7 

Q® = ZiZ®-A®/7 


} 


(B14) 

(B15) 


Similarly, 


7/u)'=i 






AoZ®' 
n , Do 


■'-‘''••(^"■'+37+1) 


Y (B22) 


The remaining quantities follow from equations (B14) and 
(B15). 


PROBLEM B 


The solution of problem B is restricted to angles ^21 for 

which — m>— (1 — r), and involves a reflected sound wave 
(X 2 

and vorticity wave behind the shock. 'Phe momentum equa- 
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tions (7) provide 


F=—%aK 


H=- 


/F' 

CZ-2 

Tv 


0K 


(B23) 


and the continuity equation (8) provides that 

(B24) 

All perturbation quantities vanish ahead of the shock, and 
shock relations (6a), (6b), (6c), and (6d) yield, respectively, 


G=L(1-Bi)-F-Ai 

Q=C,L-Ai-Kh 

K=D,L-A^ 

j^^ m{l-r)Fa2lV 

By solving the foregoing set of equations for L, 
^ 7+1 / j I 7—1 X 


(B25) 

(B26) 

(B27) 

(B28) 


L' = 


where 


+ 3 27 






1 + 


oTV 


m (1 — r) 


1 — r 


X = -((cr + m). 


1+^ r 


(B29) 
(B30) 


The remaining quantities follow from equations (B23) 
through (B27). 

PROBLEM C 

The analysis of this problem parallels that of problem A. 
1 . — The equations of motion (7) 

and (8) provide 

F=-%Kcl 
7 V 


I-- 


CI 2 

'Tv 

G 


Kp 


7l(l— /•) 

The shock equations (5) give 

G=L-F-Bx{L-Ai) 


q=C,{L-A,)-Kly 

K=D,{L-A,) 


L=.^{HAI-A2) 

rn 


These equations may be solved to yield 

L ,r, 7+1 r(l-7Kr) 


Z'=- 


=l 


V Ai^aa2^ 

where, in this case. 


1 7+1 _ 

4 7 


+ 




(7 =n (1 — r) 




(B31) 


2. The analysis is essentially identical to the 

corresponding part of problem A, except that the constants 
B 2 , Gi, and do not appear in the present case, and the 
Mach number M does not appear explicitly. Except for 
these differences, equations (Bl4) through (Bl7) may be 
carried over to the present case. The solutions for and 
may be written as 






A'lA/2 + ^fi'3 


■\ Ai A Ai^ 


i«)/ = 




■yj Ai^ A A-, 




h\ — lh/2 

WAW 


The definitions of the h’s follow those of problem A (eqs. 
(Bl9) and (B21)), with the exceptions that 


, 1— r 1 /Z \ 

ft] = I m 

r or \<7 J 


and, in this case, o- = n(l— r). The solutions for the F's and 
H’s are: 

//(« 


Fw> = . 




T^Taa^ 

iT(2) 

VAi^+/12^ 


= ^4i'''' + ^5(i“’'-0 




S- — A4(i»)'-Z) + aA,i‘^' 




V+i'+.A2‘ <7 

iJ(2) 




V +i^+ Ai^ O' 

The remaining quantities follow as before. 
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